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THE LOCAL COUNTING FUNCTION OF OPERATORS OF 
DIRAC AND LAPLACE TYPE 

LIANGPAN LI AND ALEXANDER STROHMAIER 


Abstract. Let P be a non-negative self-adjoint Laplace type operator acting 
on sections of a hermitian vector bundle over a closed Riemannian manifold. In 
this paper we review the close relations between various P-related coefficients 
such as the mollified spectral counting coefficients, the heat trace coefficients, the 
resolvent trace coefficients, the residues of the spectral zeta function as well as 
certain Wodzicki residues. We then use the Wodzicki residue to obtain results 
about the local counting function of operators of Dirac and Laplace type. In 
particular, we express the second term of the mollified spectral counting function 
of Dirac type operators in terms of geometric quantities and characterize those 
Dirac type operators for which this coefficient vanishes. 


1. Introduction 

Let M be a closed Riemannian manifold of dimension d and with metric g. Let 
E be a smooth complex hermitian vector bundle over M. As usual we denote by 
E) the space of smooth sections of E, by L‘^{M;E) the Hilbert space of 
square integrable sections equipped with the natural inner product dehned by the 
hermitian structure on the hbres and the metric measure fig on M. 

A second order partial differential operator P : C°°{M]E) — E) is said 
to be of Laplace type if its principal symbol ap is of the form ap{^) = gx{,^,0^dE^ 
for all covectors ^ G T*M. In local coordinates this means that P is of the form 

(1.1) P = -g'‘\x)didj + a’"{x)dk + b{x), 

where a^, b are smooth matrix-valued functions, and we have used Einstein’s sum 
convention. Given a Laplace type operator P, it is well known that there exist a 
unique connection V on E and a unique bundle endomorphism V G End(E)) 

such that P = V*V -|- V. Then P is called a generalized Laplace operator if U = 0. 
We will assume here that P is self-adjoint and non-negative. Thus there exists an 
orthonormal basis for E) consisting of smooth eigensections such that 

Pijij = X‘j4>j, where Xj are chosen non-negative and correspond to the eigenvalues of 
the operator P^^^. 
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A first order partial differential operator D : E) —)■ E) is said 

to be of Dirac type if its square is of Laplace type. This means that the principal 
symbol an of D satishes the Clifford algebra relations 


(1.2) ii)idi?, (^, 7] G T*M). 


Given a Dirac type operator D, we denote by 7 the action of the Clifford algebra 
bundle on E generated by the principal symbol of D, and suppose V is a connection 
on E that is compatible with the Clifford action 7 (see Section 3 for details). We 
call the triple {E, 7 , V) a Dirac bundle, and can express D uniquely as D = 7 V + 
where ^|J G C°^{M]End{E)) is called the potential of D associated with the Dirac 
bundle {E, 7 , V). In particular, D is called the generalized Dirac operator associated 
with the Dirac bundle (T^, 7 , V) if '0 = 0. We also assume D is self-adjoint, which 
means there exists a discrete spectral resolution where is 

an orthonormal basis for L‘^{M]E), and D(j)j = fijCpj for each j. Obviously, (j)j will 
be eigensections of P = with eigenvalues /i^. Therefore, using the notation from 
before Xj = |/ijj. 

Given a classical (polyhomogeneous) pseudodifferential operator A of order m G M 
acting on sections of P, the microlocalized spectral counting function iV^(/i) of D 
is dehned as 


(1.3) 


NA{li) 


jT.oiAiAMi' it 1" > 0. 

itn<0. 


Thus, is a piecewise constant function on M such that 




where 5^^. denotes the delta function on M centered at /i^. In case A is the operator 
of multiplication by a function f[x) G C°°{M) then 


where 


AC(/n) = / f{x)N^{^)d^g{x), 
J M 


(1.4) 


NM 


if/i>0 , 

lE;.<^,<oll<^i(^)llk if/i<0 


is the so-called local counting function of D. 

Let X G ^(R) be a Schwartz function such that the Fourier transform >^7 of y is 1 
near the origin and supp(t^x) C (—5, <5), where 5 is a positive constant smaller than 
half the radius of injectivity of M. It is well known (see e.g. [10, 30, 31, 42, 43, 52] 
for various special cases) that 


(1,5) (x * NMia) = B)!!"”-*-* 

j=l k=0 


{fi — 00 ). 
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This can be derived from studying the Fourier integral operator representation of 
^ Sign(D)+idi; ^-it\D\ stationary phase method. The mollihed spectral counting 

coefficients do not depend on the choice of x, and are locally computable 

in terms of the total symbols of and \D\. Note also the corresponding 

expansion for fi ^ —oo can be easily obtained from replacing D with —D: 

OO 

( 1 . 6 ) 

k=0 

Therefore, the function y * iV^ contains all the information about ±T*)}^o- 

One of the purposes of the paper is to show how to explicitly determine sz/k{A, D). 
In particular, for any bundle endomorphism F of E, we can express sz/i(F,D) in 
terms of geometric quantities such as g,'y,'Vy'lp, F. To compare, Sandoval ([43]) 
obtained an explicit expression of 32/1 (He, T*), while Branson and Gilkey ([3]) can 
also do so for £/i{F, D) whenever F is of the form /He where / is a smooth function 
on M. In the case of more general first order systems with F = He a formula was 
also more recently obtained by Chervova, Downes and Vassiliev ([7]). In the case 
of operators of Dirac-type our paper also explains the relation between the results 
of [7] and the subsequent [8] on one hand and known heat-trace invariants on the 
other. 

The mollihed spectral counting coefficients £^k{F,D) closely relate to the local 
counting function of D. Recall D is a self-adjoint Dirac type operator with spectral 
resolution {4>j, For each j we denote by = Tj(a;,|/) the Schwartz kernel 

of the orthogonal projection onto the space spanned by 0^, that is, 

= / ^j{x,y)cl){y)diig{y), 

J M 

where / e C°°{M]E), '^j{x,y) maps 0(i/) G Ey to {(j){y),(j)j{y))Ey(j)j{x) E E^. We 
then denote by = ‘hj(x) = Tj/x, x) G G°°(M; End(F')) the diagonal restriction 
of Tj. Similar to (1.5) there exists (see e.g. [30]) an asymptotic expansion 

OO OO 

(1.7) - ft) ~ E (!< ^ «i). 

j=l k=0 

where Jfk(F) G C°°(M;End(F)). This immediately implies that 

( 1 . 8 ) ^k(F,F>)= [ TiEiF^kiD)) 

Jm 

for all non-negative integers k. On the other hand, it is possible to recover ^k{D) 
from ^k{F, D) whenever explicit expressions of ^k{F, D) are known for all F. Note 
at each x G M, TrE,,,(‘hj((r)) = ||0j(a:))||E^- Thus 

OO 

(1.9) (x*Y)(l‘)~E#‘“^‘TrE.(-^‘T)W) 

k=0 

whose asymptotic expansion coefficients TrE,,,(=Sffc(D)(x)) are determined by Afk{D). 
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The mollified spectral counting coefficients s^k{F, D) also closely relate to the 
small-time asymptotic expansion of the Schwartz kernels of and De~^^ . Let 
Q : C°^{M-,E) —>■ C°^{M]E) be a partial differential operator of order m. It is 
well known (see e.g. [14, 18, 48]) that is a smoothing operator with smooth 

kernel K(t, x, y, Q, D"^) in the sense of 


(1.10) {Qe ^°^(p){x) = [ K{t,x,y,Q,D^)(p{y)dyg{y), 

J M 

where iL(t, x, ?/, Q, D"^) maps Ey to E^, and the diagonal values of K admit a uniform 
small-time asymptotic expansion 

OO 

(1.11) K{t,x,x,Q,D'^) ^ J^k{Q,D'^){x) (t-)-0+), 

fc =0 


where D"^) £ End(i?)) vanishes if k + m is odd. It is not hard to show 

for all non-negative integers k that 


( 1 . 12 ) 


•^k{D) 


T{^) ^ r(^^±i^) ’ 


where T denotes the Gamma function. This immediately implies that one is able to 
recover all the J^(Z1,Z1^) if the dimension d of M is odd and if 

d is even from all the ^k{D). To the authors’ knowledge, although how to express 
say for example k < 10, is a well-studied topic ([19, 20]), it seems that 
no explicit formula for D'^) with k odd is stated in the literature. 

Now we state the main results of the paper. Let Sub(y4) denote the principal 
and sub-principal symbol of A, respectively (See Section 4 for details), and let T^M 
denote the unit cotangent bundle of M. For any self-adjoint Dirac type operator D 
of potential ijj associated with the Dirac bundle (F, 7 , V), define ijj = 7(ej)'07(ej), 
where is a local orthonormal frame in T*M. Obviously, G End(£')) 

is independent of the choice of local orthonormal frames. Recall that d denotes the 
dimension of M. If fc = 0, then it is easy to show that 

(1-13) <{A,D) = L .a f TrE{(TA + (TA-(^D)- 

Jt*m 


Theorem 1.1. Let D he a self-adjoint Dirac type operator and let A be a classical 
pseudodifferential operator of order m on sections of E. Then 


^i{A,D) 


2 • (27r)'=' 
1 


't*m 


Trs Sub(Al) - 


2 ■ (27r)'^ Jrp* 


Tr, 


d -|- m — 1 
2 

d -|- m 


■ a A ■ Sub(D^) )-h 


M 


^Sub(y4D)- - — ■ aA- (^D ■ Sub(D^) j. 


The subprincipal symbols appearing in the above expressions depend on a choice 
of local bundle frame. The integrals are however invariantly defined as integration 
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over the co-sphere at each point yields a well dehned density on M (see Section 4). 
Note that 

(1.14) Sub( 74 i?) = Sub(A) ■ as + cta- Sub(i?) -|- ctb}, 

where B is any other classical pseudo differential operator on sections of E, {a^, <tb} 
denotes the Possion bracket between a a and as- This means that M(^, D) depends 
only on the principal and sub-principal symbols of A and D. In general, our method 
of proof can also provide an explicit formula for s^/k{A,D) in terms of the local 
symbols of A and D for each k >2. In the case when A is a bundle endomorphism 
this specializes to the following theorem. 


Theorem 1.2. Let D he a self-adjoint Dirac type operator of potential ^|J associated 
with the Dirac bundle (i?, 7 , V) and let F be a smooth bundle endomorphism of E. 
Then 


(1.15) 


^i{F,D) = 


( 47 r )'='/2 . r (|) 


TrE(F 


— {d — 2)'0, 


' M 


As immediate consequences of Theorem 1.2, one obtains 


(1.16) 


jr,(D) 


1 ip — {d — 2 )-^ 

( 47 r )''/2 • r ( f ) 2 


and 

(1.17) 


jei{D,D'^) 


1 — {d — 2)'if 

{AttYD 2 


Theorem 1.3. Let D be a self-adjoint Dirac type operator. Then ^i{D) = 0 if and 
only if D is a generalized Dirac operator. 


Note that in the special case rk(£') = 2 and d = 3 formula (1.16) implies that 
pointwise =2fi(il) is proportional to the identity matrix. Hence it vanishes if and only 
if the local coefficient tr(.ifi(il)) vanishes (see also Example 5.5). Then Theorem 
1.3 directly implies the characterization theorem in [ 8 ] (see Example 5.6 for details). 
Our result can therefore be seen as a generalization of Theorem 1.2 in [ 8 ] to higher 
dimensions and a more general geometric context. 

Finally, we explain how to explicitly determine s^k{A,D). Recall that P is a 
self-adjoint Laplace type operator with spectral resolution and A is a 

classical pseudodifferential operator of order m G M. Similar to (1.5) one has 

00 CO 

(1.18) (x*lvy(A) = ^{8l,),j,,),,)x(A-Aj)~^j4(TDA"”-'=-' (A^cx.). 

j=l k=0 

where the microlocalized spectral counting function Na{)^) of P is defined as 

(1.19) NA{X)=J2(^^L^j)- 

\j<X 
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Then it is straightforward to show for all non-negative integers k that 

(1.20) D) = 

which means to determine s^k{^i D) h snffices to do so for D^) for arbitrary 

classical psendodifferential operators A. Apart from the Fourier integral operator 
representation method introduced before, there exist several other ways to recover 
the mollihed counting coefficients ^(A, P). First, the microlocalized spectral zeta 
function ({s, A, P) is dehned by 

(1.21) C{s,A,P) = (Re(s) > d + m). 

A>o 1 

It is well-known (see e.g. [10, 52]) that ({s, A, P) admits a meromorphic continuation 
to C whose only singularities are simple poles at s = d + m — k (/c = 0 , 1 , 2 ,...) with 
residues £/k{A,P). Second, the Mellin transform of 

tr(Ae"*^) - ^ {A(j)j, (pj) {t e (0, oo)) 

A,=0 

admits a meromorphic continuation C(2s, A, P)F(s) to C whose singularities can be 
completely determined from those of ({s, A, P) and F(s). After establishing suitable 
vertical decay estimate for (^(2s, A, P)F(s), one can deduce from the inverse Mellin 
transform theorem the following widely used heat expansion (see e.g. [22, 25, 26, 
37, 47]) 

OO 

(1.22) tr(Ae-*^) ~ {m{A, P)t^ + %{A, P)tHog{t) + ^k{A,P)t’^), 

k=0 

where f —)■ O’*", and the relations between the mollihed counting coefficients and some 
of the heat coefficients can be summarized as follows: 

Case 1: If the order m of A is an integer, then 

p/ d-\-m — k \ 

• ^k{A, P) = -y- -^kiA, P) {d + m — k is positive or negative but odd); 

• ‘^fc(A, P) = 0 (d -|- m -|- 2A: < 0); 

• ‘^fc(A, P) = ^ 2 .fc!— ■ '^d+m+ 2 k{A, P) {d + m + 2k > 0). 

Case 2: If the order m of A is not an integer, then for all non-negative integers k: 

•n/ d+m —fc \ 

. i^fc(A,P) = ^^^1..^4(A,P); 

• ‘^fc(A,P) = 0. 

Thus the heat expansion (1.22) contains all the information about {. 24 (A, P)}^g. 
In exactly the same way, the following resolvent trace expansion (see e.g. [23, 25, 47]) 

OO 

tr(A(l + tP)-f) ~ ^ ^(71, P)t^^ + P)t^ log(f) + ^(A, P)t^) 

k=0 

also contains all the information about {.i 24 (A, P)}^q, where N is any complex 
number such that Re(A^) > max{d -|- m, 0}. To summarize, there exist at least four 
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different ways (mollified counting functions, spectral zeta functions, heat expansions, 
resolvent trace expansions) to determine the mollihed counting coefficients P). 

For example, using parametrix constructions in any of these methods results in the 
well known leading term 

(1.23) = TtM. 

[ZTT) Jt*M 

In this paper we will use the Wodzicki residue (see e.g [27, 34, 44, 47, 50, 51]) to 
study P). On the algebra \1/(M; E) of all classical pseudo-differential operators 

on C°°{M\E), there exists a trace called Wodzicki’s residue or non-commutative 
residue which is dehned by 

(1.24) res(A) = / ieSx{A)dx, 

J M 

where 

(1.25) TeSx{A)dx ^ [ TT{a-d{A){x,^))dS{^))dx 

is independent of the choice of local coordinates and thus is a global density on 
M, a_d{A){x,^) denotes the homogeneous part of degree —d of the total symbol of 
A, dS{^) denotes the sphere measure on If M is connected any trace r on 

\h(M; E) is a multiple of res. It is well-known (see e.g. [1, 16, 22, 33, 34, 44, 47, 52]) 
that Wodzicki’s residues are closely related to the mollihed counting coefficients as 
well as the heat coefficients. The connecting formula is ([27, 51], see also [22, (0.2)], 
[36, (1.2)], [44, (1.16)1) 

(1.26) ies{A) =-2%{A,P), 

which means 'ifo{A,P) is independent of the choice of P. As for other coefficients, 
it is known for all non-negative integers k that (see e.g. [16, Thm. 5.2]) 

/I \fc-i-i 

(1.27) «i(A,P) = YJ^.res(2lF'=), 
and (see e.g. [34, Prop. 4.2], [47, P. 106]) 

(1.28) Mi(A, P) = res(AP^^). 

Thus combining (1.20) with (1.28), one gets 

(1.29) j4(A, D) = res(A ^^*^* |P|^-'^-"^-^), 

which is the main tool for us to establish Theorem 1.1. We should mention that any 
of the previous four methods can also be used to derive Theorem 1.1, but it seems 
that these methods do not provide such a clear interpretation of £fk{A, P) for k >2. 

This paper is structured as follows. 

Section 2 hrst reviews the construction of the Fourier integral operator 
and the stationary phase expansion of then provide as consequences complete 

proofs of the singularity structures of the spectral zeta functions, the heat expansions 
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as well as the resolvent trace expansions. As applications, (1.27) and (1.28) follow 
from snitably applying formnla (1.26). The local connting fnnction of operators of 
Dirac and Laplace type is also stndied in the last part of this section. 

Section 3 hrst reviews the concept of the Dirac bnndle and the associated Dirac 
type operators, then gives a direct proof of Theorem 1.2 by appealing to the Bochner- 
Weitzenbock techniqne in Riemannian geometry. 

Section 4 proves Theorem 1.1 and Theorem 1.2 by applying the Wodzicki residne 
method. The connections between certain Wodzicki residnes and the snb-principal 
symbols of classical psendodifferential operators are also discnssed. 

Section 5 characterizes Dirac type operators with vanishing second expansion 
coefficient. As an application, the so-called massless Dirac operators are stndied. 

Thronghont the paper we assnme on smooth sections of E, 

• F is a self-adjoint Laplace type operator with spectral resolntion {(j)j, 

• F is a self-adjoint Dirac type operator with spectral resolntion 
nnless otherwise stated in the last section, 

• A is a classical (polyhomogeneons) psendodifferential operator of order m, 

• Q is a. partial differential operator of order m, 

• F is a smooth endomorphism. 


2. Mollified counting coefficients I: qualitative theory 

This section hrst reviews the constrnction of the Fonrier integral operator 
and the stationary phase expansion of y * N'^, then proves many conseqnences. 


2.1. FIO method. Formnla (1.18) essentially is Proposition 2.1 in [10], Corollary 
2.2 in [30], Theorem 2.2 in [43] and Proposition 1.1 in [52], except they either consider 
scalar operators or assnme A is of order zero. Recall that y G A^(M) is chosen so 
that = 1 near the origin and snpp(t^x) C (—5, 5), where 6 is smaller than half 
the radins of injectivity of M. If t is snfhciently small, say |f| < (5i < 5, then locally 
the integral kernel (Ae“‘*'^^^^)(f, x, y) of the operator is well known to have 

the form 


(Ae '*-^'^')(f,x,|/) = 


(2vr) 




where a is a classical (matrix-valned) symbol of order m. The scalar-valned phase 
fnnction 9(t,x,y,^) = K{x,y,^) — tapi/ 2 {y,^), where K{y,y,^) = 0, was introdnced 
by Hormander ([29]). It is also known that tr(Ae“'*'^^^^) is smooth in (—5, 5)\{0}, 
so we introdnce a cnt-off fnnction g G ^^'(M) satisfying g{t) = 1 if |f| < y and 
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supp(^)) C (—(5i,(5i). Using integration by parts one gets 


{x*N'^)iX) =— / (^x)(t)(p(t) + 1 - p(t))tr(Ale 

./h 


iAi 




( 2 ir: 


d+1 



'M JR<i 


{’^x){t)Q{t)TT{a{t, y, y, ^))e 


+ o(A-°°) (A^cx)), 


where o(A °°) is short for o(A for all positive integers h. Now proceed as in [10] 
via the stationary phase method to get (1.18). 


2.2. Finite heat expansions. It is easy to prove that 

(2.1) (T, yPi) = {p*T, + o{n {t ^ 0+), 

where T is a tempered distribntion on M, (p G =.5^(M), ptW = and p G J^{R) is 

chosen so that ^p = 1 near the origin. Actually, this formula appears in equivalent 
forms in [ 10 , 28, 52], so its proof is omitted. As an application, one easily gets 

( 2 . 2 ) {N'^, AV(tA)) = (x * N'^, AV(tA)) + o(f-) {t ^ 0 +), 

where h is an arbitrary non-negative integer. Note the left hand side of (2.2) just is 
tr(AP 2 (p(t\/P)). Now we claim the right hand side of (2.2) is of the following form 

poo 

(2.3) j4(A, P) ■ / A'^+™+'^-^-V(A)dA ■ + o(tM-™-i). 

k<d-\-m-\-h ^ 

To this end we first decompose 
(X' » K, AV(n)) = / (x» K)(X) ■ AV((A)<iA + 

J — oo 

/ CX) 

M:(A,P)A'^+”^-^-^ • AV(fA)dA + 

K<^a-trri-tn 

{ix * K)iX) - Y P)A'^+—fc-i) . AV(fA)dA 

k<d-\-m+h 

= ai{t) + a2{t) + a3{t). 

Since supp(A^j^) C [0, oo),it is easy to verify that (x*Ad^)(A) = o(|A|“°°) as A —?■ — cx). 
This implies Q!i(f) =0(1) as f —)■ O'*". Note also 

02 (t) = Y ■^k{A,P)- 

k<d-\-m-\-h 

For simplicity we introduce 

/(A) = ((x»/V^)(A)- i4(4DA‘‘+”-‘-')-A'‘ (A>0) 

k<d-\-m-\-h 

and note a^it) = f{X)(p(tX)d\. To prove the claim we have two cases to consider. 

Case 1: Suppose d -|- m -|- h > 0. Let k be the unique non-negative integer such that 


k d-\-m-\-h—k—l 


p{X)dX ■ t 


k—d—m—h 
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k<d + m + h<k + l. Let I3 = d + m + h — k — 1 G (—1, 0], which implies there 
exists a constant Ci such that for all A G (0,1), |/(A)| < CiA^. According to (1.18), 
there exists another constant C 2 such that for all A > 1, |/(A)| < C 2 A^“^ < C' 2 A“^. 
So as f —)■ 0+, 

/ I noo 

\^\ip{t\)\d\ + C 2 j \-^\ip{t\)\d\ 

pt poo 

= Ci A^|(p(A)|dA-r (^+^)+^2 / A-^|<p(A)|dA 

Jo Jt 

= o(r''’+«) + 0(|logMI) 


Since —(/3 + l) = k — d — m — h = \rri \ — m — 1 < 0, we are done in the hrst case. 
Case 2: Suppose d + m + h < 0. Obviously, there is a constant C 3 such that for all 
A G (0,1), |/(A)| < C 3 . According to (1.18), {x * N'^){X) = 0(A'^+"*“^) as A —)■ cxo. 
Thus there is a constant C 4 such that for all A > 1, |/(A)| < C^X-\ So as f ^ 0 +, 

pi poo 

Mt)\<C3l \^itx)\dx + cj x-^\^itx)\dx 

= 0(l) + 0(|log(f)|) 

= C»(|log(t)|). 


But note [m] — m — 1 < 0, we are also done in the second case. 
To summarize, we have shown for all G ^(M) as t —)■ O’*' that. 


(2.4) tr{APiip{t^)) 


/ . ^d^m+h-k ) 

k<d-\-m+h 


2 . 3 . Mellin transforms. The Mellin transform of a continuous function / on 
( 0 , 00 ) is the function (M/)(s) of the complex variable s, given by 

poo 

{Mf){s) = / f{t)f-^dt 
Jo 

whenever the integral is well-dehned. An open strip n(/ 3 i,/ 92 ) = {s G C : /3i < 
Re(s) < 132 } is called a basic strip of M/ if the integral is absolutely convergent in 
that strip. Given meromorphic functions u, v dehned over an open subset If of C, 
we denote m x n (s G If) to mean m — n is analytic in If. For example, assume / is a 
continuous function on ( 0 , 00 ) satisfying f{t) = o(t~°°) as f —?■ 00 and assume there 
exist real numbers wq < Ci;i < • • • < wat such that 


N-l 

/(t) = + 0 (f‘"^) (t^ 0 +). 

k=0 


(2.5) 
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Then one may easily check that n(— cjq, cxo ) is a basic strip of Mf, and Mf admits 
a meromorphic continuation to n(— oo) such that 


N-l 


( 2 , 6 ) 


(M/)(s) X 




k=0 


S + Uk 


(s G n(— cjat , cxo )). 


Actually, for any s G 11 ( 71 , 72 ) with —ojn < 71 < 72 < cxd , one has 


N-l 


{Mf){s) = Y^ ^ + I f{t)f ^dt, 


N-l 


k=0 


s U!k Jo 


k=0 


which immediately implies ( 2 . 6 ) as well as an upper bound estimate: 

(2.7) (M/)(s) = 0(l) (s G n(7i,72), |s| ^ 00). 

Lemma 2.1 ([15]). Let f be a continuous function over (0, 00 ). Assume there exist 
real numbers (Ji < (32 < ^3 < (3a such that 


• Il{(3z, (3a) is a basic strip of Mf, 

• Mf admits a meromorphic continuation to n(/ 3 i,/ 34 ) with finite poles there, 

• M/X (,1^ (.Gn(/32,/34)), 

• Mf is analytic on Re(s) = (32, 

• {Mf){s) = 0(|s|-2) (s e n^i,/34, |s| ^ 00 ). 

Then 

f(t) = c „,,( hA (-( iog ( y ) + o ( r *) ((^ o +). 


2.4. Spectral zeta functions. In this part we deduce the well-known singularity 
structures of the spectral zeta functions. Let h be a positive integer such that 
d -|- m -|- h > 1 , and dehne 

( 2 . 8 ) fh{t) = tr(APte-*^) = ^ {Acfj, {t > 0 ). 

Aj>0 

Now we have two ways to study the Mellin transform of fh- First, we list {Aj}^i in 
non-decreasing order and thus by Weyl’s law (see e.g. [2]), limj^oo exists and is 
positive. It is easy to see (see e.g. [4, 48]) that there exists a constant C = C(A, P) 
such that \{A(f)j,(j)j)\ < CA™ whenever Aj > 0. Consequently, 

• fh is smooth over (0, 00 ) with fh(t) = o{t~°°) as t —)■ 00 ; 

• the spectral zeta function ({s, A, P) is analytic in n(d -|- m, 00 ); 

• Mfh has a basic strip n(^Ti|±k^ qq) in which [Mfh){s) = ({2s —h, A, P)T{s). 


Second, according to (2.4) fh is easily seen to have the form (2.5): 


(2.9) 


fh{t) = 


N—l d+m+h—k ' 


k — d — m — h 


{A,P)-t —^— +0{P^) (f^0+). 


A:=0 


2 
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where AT = Nh = d+h+\m\-l > 1 , wjv = < 0 , 0 ;^ = a;jv-^ = 

(/c = 0,1,..., N — 1). So by (2.6), Mfh initially analytically defined in oo), 

admits a meromorphic continuation to n(— oo) in which 


( 2 . 10 ) 


(MA)(s) 


A -1 


E 


r( 


d-\-ra+h—k 

2 


) 


2 


MAP) 

_ dAm+h—k * 

2 


Considering in n(— wat, oo) the Gamma function F is analytic and has no zeros, it 
is easy to see that ({s, A, P) initially analytically defined in n(d + m, oo), admits a 
meromorphic continuation to n(— 2 ci;Ar — h, oo) in which 


( 2 . 11 ) 


Q(s,A,P) 


(MA)(^) „dC* MA,P) 


Letting h —)■ oo we get that ({s,AiP) admits a meromorphic continuation to C 
whose only singularities are simple poles at s = d + m — A; (A; = 0,1, 2,...) with 
residues P). 


Remark 2.2. Note that if P is a smoothing operator on sections of E, then 
(’(s,P,P) is an entire function on C, and consequently, s^k{A,P) = ^k{A + B,P) 
for all k. 


For any g G M, we let P^ : E) —)> E) be the operator defined by 

functional calculus of P if g is non-negative and by sending 0 to X]aj>o 4’j)4’j 

if g is negative. By a classical result by Seeley ([48]) the operator P"^ is a classical 
pseudodifferential operator of order 2 g. 

Proposition 2.3. For any real number q, s^k{AjP) = ■^k{AP^,P) holds for all 
non-negative integers k. 

Proof. Note AP'^ is a classical pseudo differential operator of order m + 2g. Thus 
C(s,y4P^,P) admits a meromorphic continuation to C whose only singularities are 
simple poles at s = d -|- m -|- 2 g — A; (A; = 0, 1 , 2 ,...) with residues a 24(^P'^, P)- But 
noting that C(s, AP"^, P) = M — 2g, A, P), we must have s^k{AP'^, P) = .(^k{A, P). 
This finishes the proof of the proposition. □ 

Remark 2.4. To determine s^k{A, P) we can assume without loss of generality that 
A is of order zero as s^k{AP~^, P), which is s^k{A, P) by Proposition 2.3. 

2.5. Full heat expansions. In this part we provide another approach to the heat 
expansion (1.22). Define a smooth function / over (0, cxd) by 

(2.12) f{t) = ^ {A(j)j, f)j)eM = tr(24e"‘^) - ^ 0^) (f > 0). 

A j O 0 Xj=0 

Obviously, n(max{^^, 0 }, cxo) is a basic strip of Mf which admits a meromorphic 
continuation to C with {Mf){s) = ({2s, A, P)r(s). Recall <(’(•, A, P) is meromorphic 
on C whose only singularities are simple poles at s = d-|-m — A; (A; = 0,1, 2,...) with 
residues £^k{A, P), while T is meromorphic on C whose only singularities are simple 













DIRAC AND LAPLACE TYPE OPERATORS 


13 


poles at s = —k {k = 0,1, 2,...) with residues ^ . Hence to establish (1.22) via 

Lemma 2.1 it suffices to prove 

(2.13) {Mf){s) = 0{\s\~‘^) (sGn(—n,n), |s| —)■ oo) 

for all positive integers n. To this end we hrst assume s G n(— n,n) {n G N), then 
let h > n + 2 be a large enough integer so that one can use ( 2 . 11 ) to get 

(2.14) (MS)(s) = ({2s,A,P)T{,8) = 

' ' i=0 

Finally by considering s + h G n(2, n + h), (2.13) follows immediately from suitably 
applying (2.7) to uniformly bound {Mf 2 h){s + h) above for all s G n(— n,n) with 
|s| —?■ oo. This hnishes the proof of (1.22). 

The precise relations between the mollihed counting coefficients and some of the 
heat coefficients as mentioned in the introductory part are left to the interested 
readers to verify. For completeness we record these relations as a proposition for 
later use. We also remark that the above proof actually yields 

(2.15) C(2s, A, P)F(s) = 0(|sr°°) (s G n(-n, n), |s| ^ oo) 
for all positive integers n. 

Proposition 2.5. If the order m of A is an integer, then 

p/ d-\-m — k N 

• IiSk{A, P) = -^— • s^k{A, P) (d + m — k is positive or negative but odd); 

• ^k{A,P) = 0 (d + m + 2k<0); 

• ^k{A, P) = ^ 2.^— ■ •^d+m+2k{A, P) (d + TTl + 2k > Oj. 

If the order m of A is not an integer, then for all non-negative integers k: 

1-'/ d-\-m — k \ 

. ^,(h,p) = ^^^-34(AP); 

• %{A,P) = 0. 

Remark 2.6. There exists an even stronger estimate than (2.15) in [26]. But to 
apply the inverse Mellin transform theorem, a rough estimate like (2.15) is sufficient. 


2.6. Resolvent trace expansions. In this part we study the short time asymptotic 
expansion of the resolvent trace {f^^\t) =) tr(H(l +tP)“^P), where iV G C is such 
that Re(A^) > max{d + m, 0}. Let denote the Mellin transform of 


tr(H(l + tP) ^/2) - {Af)j, (fj) (t > 0). 


Aw=0 


It is easy to verify that has a non-empty basic strip n(max{^^, 0 }, j 

which = C(2s, A, P)B{s, y — s), where we recall the Beta function 

F(a)F(/3) 


m 


B{a,(3)= / (1 — t)" ^t^ ^dt = 


u 


OL — 1 


(1 +n) 


a+g 


du = 


F(a + /3) 
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is defined for Re(Q!) > 0 and Re(/5) > 0. Hence admits a meromorphic 

continuation to n(— oo, which 


(2.16) 


mW(s) = C(2s, 4 


For any real /3 in the strip n(max{^^ 
easy to see that 

sup |r(— -s)|< 

s£n(—n,/3) 


, 0}, any positive integer n, it is 


max r(t) < oo. 


Consequently by (2.13), 

(2.17) = 0{\s\-'^) (s G n(-n,/3), |s| ^ oo), 


which immediately implies 

OO 

(2.18) /(^)(t) ~ )(AP)t"log(t) + 

k=0 

via Lemma 2.1 as the singularities of ({■,A,P) and F are known to be completely 
determined, so is The precise relations between the mollihed counting coef- 

hcients and some of the resolvent trace coefficients are summarized in the following 
proposition: 


Proposition 2.7. If the order m of A is an integer, then 




(A,P) = 


d-\-m — k j N-\-k — d—7n j 


r(f) 


J 2 

negative but odd); 

^1^\A, P) = 0 Cd + m + 2fc < o;,- 


• £^k{.A, P) (d + m — k is positive or 


. )(H, P) = • ^d+m+2k{A, P) (d + m + 2k> 0). 

If the order m of A is not an integer, then for all non-negative integers k: 

/ \T\ -p/d+m —p/iV+fc —d —m \ 

. P) = . S=f=^ ■ i^k(A, P); 


r<^\A,P) = 0. 


2.1. Wodzicki residues. This part is devoted to the proofs of (1.27) and (1.28) 
and their consequences. Both formulae are known and we state them here as two 
propositions for completeness. Also, we refer the readers to the original references 
for the proof of the crucial formula (1.26). 

Proposition 2.8. s^k{A,P) = res(AP*' 2 "“), 

Proof. Since AP’° 2 "" is of integer order k — d and d + (/c — d) + 2 • 0 = /c > 0, it 
follows hrst from Proposition 2.5 then from Proposition 2.3 that 

^ . k — d — m , . , . k — d— m , , , , , 

-2%{AP^,P) = s^k{AP^,P) = s^k{A,P). 
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But according to (1.26) we also have 

^ . _ k — d — m _, , , _ k — d — m , 

—2^q{AP 2 ,P) = Yes{AP 2 ). 

Combining the above two formulae yields the desired result. □ 

Proposition 2.9. ^k{A,P) = •res(74P^). 

Proof. If A is of integer order m > —d — 2k, then it follows hrst from Proposition 
2.5 then from Proposition 2.8 that 

Vt{A, P) = P) = A-2— . ies{AP‘‘). 

If A is not of integer order or A is of integer order m < —d — 2k, then it follows from 
dehnition that res(y4P^) = 0 and from Proposition 2.5 that ^k{A,P) = 0. This 
hnishes the proof. □ 

Next we study the property £/k{A, P) = 0 for certain operators A G 
We denote by the space of odd-class pseudodifferential operators and 

by E) the space of even-class pseudodifferential operators on sections of E, 

that is, A G if in local coordinates its symbol J2j>o^m-j{A) {m G 1) 

satishes 

(2.19) a^-j{A){x,-f) = {-l)^~^am-j{A){x,f) 
for all X, f and j, while A G Tg,^g^(M; E) if in local coordinates 

(2.20) am-jiA){x,-f) = (-l)'”"^+Vm_j(Al)(x,0 

for all x,f and j. For example, any partial differential operator is odd-class while 
P^/^ is even-class. An operator A G E) U E) is said of regular 

parity class if its parity class agrees with that of d. It is easy to verify that res(A) = 0 
if A is of regular parity class. 

Proposition 2.10. Let A G P) U P). If the parity class of A 

agrees with that of k — m, then s^k{A, P) = 0. 

Proof. We dehne a map r by sending odd-class pseudodifferential operators to 0 and 
even-class ones to 1. Then (see [41, Prop. 1.11]) t{AP 2 ^) = t{A) + k — d — m 
(mod 2). Thus if the parity class of A agrees with that of A; — m, or equivalently, the 
parity class of t{AP 2 ) does not agree with that of d, or further equivalently, 

the parity class of AP 2 agrees with that of d, then AP 2 is of regular parity 
class and consequently by Proposition 2.8, M:(A, P) = res(AP 2 ) = 0. □ 

Remark 2.11. The heat expansion (1.22) can be written in a more rigorous way as 

00 

(2.21) tr(Ae-*^) ~ 5^4(A,P)t"^ + (^^(A, P)t'^ log(f) + ^k{A,P)t^), 

j'eA fc=o 

where A denotes the set of non-negative integers j such that is not a non¬ 

negative integer. In this notation system A§j{A,P), ^k{A,P) and ^k{A,P) are 
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uniquely determined. According to Propositions 2.5, 2.8, 2.9, ^j{A, P) and P) 
are certain Wodzicki residues. In most situations, ^k{A,P) are not Wodzicki 
residues but operator traces. For the simplest example that A is smoothing op¬ 
erator, one can (regard A as of any real order and thus) deduce from (2.21) that 

CO 

(2.22) tr(Ae-‘^) ~ ^ ^k{A, P)tK 

k=0 

In this situation it is so natural to expect tr(A) = ^o{A,P). Actually, if any of the 
following three cases happens, then ^o(A, P) just is the Kontsevich-Vishik trace of 
A which is independent of the choice of P (see e.g. [24, 34, 38]): 

• m < —d] 

• m is not an integer; 

• A is of regular parity class. 

Since the Kontsevich-Vishik trace agrees with the standard L^-operator trace on 
trace class operators, the previous expectation is conhrmed. Furthermore, ii A = Q 
is a partial differential operator and if d is even, then one can express ^o{A, P) as the 
Wodzicki residue of certain pseudodifferential operator with log-polyhomogeneous 
symbol (see [39, 40, 46] for details). 

2.8. Dirac operators. This part is devoted to the proofs of (1.5), (1.20) and (1.29). 
First, we introduce a classical pseudodifferential operator B of order m by dehning 
B = ^ Sign(D)+idj; ^ gjj^ce supp(A'^ — N'^) C (—cxo, 0], we get (see also Subsection 2.2) 

(X * -{X* = o(/i"°°) oo). 

By appealing to (1.18) with A,P,X,Xj replaced respectively by B, we 

obtain as /r —)■ cxo that 

CX) 

(X * ^b)(/^) = - /^i) 

3 =^ 

oo 

OO 

k=0 

which proves (1.5) as well as £^k{A, D) = s^k{B, D^) for all non-negative integers k. 
Obviously by setting A = F, (1.20) is a special case of £/k{A, D) = s^k{B, D‘^). Now 
it follows from Proposition 2.8 that 

(2.23) £/k{A, D) = res(A ^^^^^^^^ ^ ^ res(A:^-^^|D|^-‘^-™-i), 

which proves (1.29). 
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Next, we give an equivalent expression for D). By Proposition 2.3 one gets 


(2.24) 




^k{A,D^) + s^k{AD,D‘^) 
2 


Similar to the microlocalized spectral zeta function (1.21), one can introduce the 
microlocalized spectral eta function (see e.g. [3, 21]) by dehning 

(2.25) r]{s,A,D) = ^ • Sign(/ij) (Re(s)>d + m). 

But it is easy to see that 

(2.26) vis, Al, D) = Cis, Al • Sign(Zl), 

So according to Proposition 2.8, 77(5, A, D) admits a meromorphic continuation to 
C whose only singularities are simple poles ais = d + m — k [k = 0,1,2,...) with 
residues res(RZl|Zl|^“'^“”*“^). 


2.9. Local counting functions. This part is devoted to the proof of (1.12) and 
its consequences. Recall Q denotes a partial differential operator of order m acting 
on C°°iM; E). Combining (1.10) ~ (l-H) yields 


(2.27) 


tr(Qe 


-tD'^ 


k=0 


TiEi^kiQ:D^))-t- 


{t —)■ 0^), 


'M 


which compared with the heat expansion (1.22) gives 


• ^kiQiD'^) = iM'^^^i^kiQ, {d + m — k is positive or negative but odd); 

• %iQ,D^) = 0. 


According to Proposition 2.5 we also have 

• ^k{Q,D‘^) = - 2 —^M:(Q,-D^) {d + m — k is positive or negative but odd); 

. ^k{Q,D^) = ^-^^^d+m+ 2 kiQ,D^) {d + m + 2k>0). 

Thus for all non-negative integers k, 

( 2 . 28 ) f 

t [ -2- ) dM 

By combining (1.8), (2.24), (2.28) and noting J^k{FQ, D^) = FJ^k{Q,D‘^), we get 


(2.29) 


TTEiF^kiD)) = s^kiF,D) = 


IM 


= / TrE(F( 

J M ^ 


^k{F,D^) + ^k{FD,D‘^) 
2 

M’k{UE,D^) , J^k{D,D^) 


P(V) 


+ 


r( 


d-\-l—k ^ 


-)). 
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which proves (1.12) as F can be any smooth endomorphism of E. To be precise, 
one can set F to be the adjoint of 

then apply (2.29) to conclude (1.12). 

Next, we study the local counting function of Laplace type operators but hrst 
have a look at that of Dirac type operators. According to (1.7), 


$iX(/i T h) - ^ ^^k{±D) (/i ^ oo), 

j=l k=0 

which easily implies (see also [30]) 




—CX) . 


t=l 


fc =0 


As an application of ( 1 . 12 ), one gets ^k{D) + ^k{—D) = . Thus 

^ V 2 ' 


2 ^$,x(p-1/1,1) ~ 2 ^/i -- (/i^oo). 


(2.30) _ _ 

i=l k=0 ■‘■''"2 

Generally for Laplace type operators P, one can similarly establish that 


(2.31) 


OO 


CX) 


5^4.,x(A-A,)~ 5^ a 

j=l k={) 


d—k—1 


2jrfc(id^,p) 

r(d^) 


(A —)■ CX)). 


3. The Bochner-Weitzenbock technique and a direct proof of 

Theorem 1.2 

In this section we will introduce some basic notations about Dirac operators and 
give a direct proof of Theorem 1.2 by employing the Bochner-Weitzenbock technique 
in Riemannian geometry (see e.g. [2, 17, 19, 35]). 

3.1. Bochner-Weitzenbock technique. We begin by reviewing the concept of 
the Dirac bundle (P, 7 , V) (see e.g. [35]), where 7 is a Clifford module structure on 
E and V is a connection on E that is compatible with 7 . 

First, let Cl(TM) be the universal unital complex algebra bundle generated by the 
tangent bundle TM subject to the commutation relation X*Y+Y*X = —2g{X, Y), 
where X, D G C°°{M\ TM) and * denotes the algebra operation. A Clifford module 
structure on E just is a unital algebra morphism 7 : Cl(TM) —)■ End(P). Given a 
connection V on the Clifford module E = (P, 7 ), 7 V = 7 (efe)Vej. is a well-defined 
hrst order partial differential operator acting on E), where {ek}f=i denotes 

any orthonormal basis of T^M at each point x G M. Next, one can naturally extend 
the Riemannian metric and connection V on TM to Cl(TM) with the properties 
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that: the extended connection denoted also by V preserves the extended metric, 
and 

(3.1) Vx(a */d) = (Vxa) */d + a * (Vx/d) 

for all X G C°°{M;TM), a,(3 ^ Cl(TM)). We then call a connection V on 

the Clifford module {E, 7 ) compatible with 7 if 

(3.2) Vx(7(«)0) = 7(Vxa)(/> + 7(«)(Vx0) 

for all X G C°^{M;TM), a G C°°(M; Cl(TM)), 0 G C°^{M;E). It is well-known 
([3]) that on a Clifford module there always exists a compatible connection. Such a 
triple {E, 7 , V) is called a Dirac bundle. With the concept of Dirac bundle available, 
it is easy to see that any operator D of Dirac type acting on E) can always 

be written as D = 7 V -|- ^|J, where (£', 7 , V) is some Dirac bundle with 7 uniquely 
determined by the principal symbol ao of D, -0 G End(£')). We call the 

potential of D associated with the Dirac bundle (E, 7 , V). 

Let (£', 7 , V) be a Dirac bundle and let = 7 V -|- -0 be a Dirac type operator 
with potential 'ip- The Bochner-Weitzenbock technique permits us to hnd a unique 
connection V■^|} on E and a unique V^p G End(£')) such that 

where denotes the connection Laplacian generated by In a series 

of papers by Gilkey and his collaborators (see e.g. [3, 19, 20]), Vp, is always written 
as a function explicitly depending on and a few others. Our feature of proving 

Theorem 1.2 is to first express Vp^ as a function of 7 , V, ip, then treat the hrst two 

variables as dummy ones. 

Given a Dirac bundle (E', 7 , V), we dehne a connection V on End(E) by 

(3.3) {Xx^)i<P) = Xxi^i^P)) ~ ^{Xx(p) 

for all X G C°°{M;TM), oj G C'°°(M; End(E)), (p G C°°{M;E). It is easy to verify 
that V is compatible with 7 in the following sense: 

(3.4) (Vx)(a; 7 (E)) = (Vxa;)7(l^) + uj'l{XxY) 

(3.5) (Vx)(7(h^)^) = 7 (h^)(Vxu;) + 7 (VxX)u; 

where X,u},(p remain the same meanings as before, Y G C°°{M]TM). 

We also need to fix a few notations. If (x^,..., x'^) is a system of local coordinates 
on M, then let {di = denote the local coordinate frames of TM. We use the 
metric tensor to identify the tangent and cotangent bundles TM = T*M, which 
means locally dx^ = g^'^di, where is the inverse of the matrix {gij = g{di,dj)). 
The curvature tensor of any connection V on E is denoted by , that is, R^y — 
[Vx, Vy] — V[x,y] for X,Y E C°°{M;TM). Let TL denote the Christoffel symbols 
of the Riemannian connection V, that is, Xa-dj = ^ijdk- The connection Laplacian 
locally is of the form —g^^{Xo^X^. — LLVaj,)- Finally, we denote = 

U{dx ^^,..., dx*") and 17*1,..., c?*„) for any map U defined on the n-fold 
product of TM. 
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Proposition 3.1. Let D he a Dirac type operator of potential if associated with the 
Dirac bundle (i?, 7 , V). Let L = L^ : TM ^ End(£') denote the map defined by 

j{X)if + if^{X) 


(3.6) 


L{X) = 


Then = X — L is a connection on E and = D'^ — e End(E)). 

Locally we have 


(3.7) 


^7 = ViV'] + L^Li + 


Proof. Let L be naturally identified with 

C~(M; TM) X C°°(M; E) 4 C°°(M; E), 

which is C'°°(M)-linear in both variables. Thus V = V — L is a connection on E. 
According to (3.3) we have 

= -g^^iXiXj - XiX, + XiXj - V,V,) + g^^T^Xk - Xk) 


= -g^^X,L^-g^^LiX,+g^%L 


= -g^^{{X,Lfi + LjXi) - g^^LfiX, - Lfi + g^^T^L^ 
= _ 2 L*V, - g^K^iLj) + VU + g^%Lk. 

On the other hand, according to (3.2) and (3.3) we have 
D^ = {-i^X, + if){-i^Xj+if) 

= f 7'V.V, - + i^iy.if) + 74V. + if-t^Xj + if^ 

— VjVi I ri “ 


— 'y 'y'f 


+ 47^’^ + + f (V4) + (74 + V^4)V. + 

~ ~ ~ L^ _ 

= + g^^T%Xk + 74"'^ + 4(V.i/^) + (74 + V'4)V. + if^ 

= A^ + (74 + V’4)Vi + 74 ''^ + 4 (V4) + V'"- 

Consequently, combining the above calculations yields 

D^ = A^ + 7 * 7 ^'^ + ^fiXiif) +if^- g^\X,Lj) + VL, + g^^T^^L^ A 

Next let us simplify the expression of 1^. According to (3.4) and (3.5) we have 


-g^^{X,Lfi = -g 


i,-74V4) + rfj74 + (V4)7i + ^ ij'iPik 


= -4^ 


7 74V4) + (V4)7j 


7*(Vi'0) + (V.'0)7* 


- 


2 
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which implies 


2 7*(ViV') + 


_ 

y* = +y- 

This finishes the proof of Proposition 3.1. 


+ VL, 


□ 


3.2. Direct proof of Theorem 1.2. Let denote the Dirac type operator of 
potential ip^ = ip — eF associated with the Dirac bundle {E, 7 , V), where e is a real 
parameter. In particular D = Dq. Let K = be as determined by 

Proposition 3.1. Locally, 


V, = ^YYRfj + \b\ + (L*.)‘(L».), + y. 


Thus 


d 

de 


e=0 


Tr(K) = Tr((L_^)*(L^), + {L^)\L_f)P) - 2Tr(Fi/.), 


which following simplification via the facts 7 * 7 * = —d and pj = 7*'07i = li'P’Y gives 

d 


(3,8) 


de 


e=0 


Tr(K) = Ti{F{{d-2)iP-iP)). 


But it is known (see e.g. [3, Lemma 2.1], [19, Thm. 3.3.1]) that 


(3,9) = 


<:=0 JM 


Tr(rIdE-6D,), 


where r is the scalar curvature of M. By combining (3.8), (3.9) and by considering 
the fact J^i{FD,D‘^) = F, D"^), one gets 


' M 


Tt{F ■ = 


( 4 ^)(i /2 


' M 




Consequently, 

(3.10) 


^i(D,D 2 ) = 


Ip — {d — 2)ip 


( 4 ^)(i /2 2 

as F can be any smooth endomorphism of E. By (1.8) and (1.12) (see also (2.29)), 


.TT)= / Tr(F(i^ 

JM ^ t 




+ 


2 . r(f) 


-))■ 


But D'^) = 0 (see e.g. [18]), which completes the proof of Theorem 1.2. 


Remark 3.2. It is known (see e.g. [18, 48]) that (1.10) and (1.11) still hold even if 
without assuming the self-adjointness of D if the integral kernels are appropriately 
interpreted. In this general situation our direct proof of Theorem 1.2 can also yield 
(3.10) with no modification at all. 
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4. Mollified counting coefficients II 


In this section we will express {A, P) in terms of principal and sub-principal 
symbols of A and P, and prove Theorems 1.1, 1.2. We refer the readers to [7, 8, 10, 
12, 30, 31, 42, 43] for some other explicit two-term asymptotic expansions in various 
situations. 


We begin by reviewing the dehnition of the sub-principal symbol of a classical 
pseudodifferential operator. 

Recall M is a closed manifold over which the space of smooth sections of the 
half-density bundle is denoted by An operator H : —)■ 

C°°{M,Qi/ 2 ) is called a classical pseudodifferential operator of order m if on every 
coordinate patch its total symbol admits an asymptotic expansion 0 

with homogeneous of degree m — j in It is well-known that the sub-principal 
symbol (see e.g. [11, Section 5.2]) of H dehned in local coordinates by 

(4.1) Sub(/l) = .2> + 5.^. 

transforms like a homogeneous smooth function of degree m — 1 on T*M\0 under 
change of charts. 

Let A : C°°{M) —)■ C°°{M) be a classical pseudodifferential operator of order m. 
In local coordinates we denote by {j = 0,1, 2,...) the homogeneous 

part of degree m — j of the total symbol of A^. Given a Riemannian metric g on M 
one can identify C°°{M) with Q 1 / 2 ) via the map 




where is the Riemannian density (called also metric measure) on {M,g). Then 
H = <l>oAo<l>“^ is a classical pseudo differential operator of order m on half-densities. 
On local coordinates let G = det{gij). Note H = AG~^I^ can be regarded as 
the product of three (local) pseudodifferential operators with corresponding total 
symbols G(a:)^A^ G{x)~^^‘^. Thus one can deduce from the product 

rule of pseudo differential operators that 


Sub(i7) = ^ 


avf i_ 

dx^dS,k 2 


cl(log \/G) 
d^k dx^ 


on local coordinates. This means the sub-principal symbol of A locally dehned by 
([42, Remark 2.1.10]) 


(4.2) 


Sub (A) = 


2 dx^d^k 


1 cl (log VG) 

2 d^k dx^ 


is a homogeneous smooth function of degree m — 1 on T*M\0. 


^Note we also denote by am-j{A) in the previous sections. 
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Now let A : C°°(M; E) —)■ E) be a classical pseudodifferential operator of 

order m. For a fixed local bundle trivialization of E over a coordinate neighborhood 
U C M, one can naturally identify : C^{U-,E) —)■ C°°{U-,E), the restriction of 
A onto U, with a matrix of classical pseudodifferential operators A^^^ : 

C°°{M) —>■ C°°{M) of order m, where r denotes the rank of E. This implies that the 
sub-principal symbol of A dehned by (see also [32]) Sub(y4) = (Sub(y4^j,))i<^,i,<r is a 
homogeneous smooth Mat(r, C)-valued function of degree m — 1 on T*U\0. On local 
coordinates Sub(y4) is still of the form (4.2). We should remember, however, that the 
dehnition of the sub-principal symbol depends on the choice of local frames for E and 
Riemannian metric on M. If the principal symbol is a multiple of the identity then 
the sub-principal symbol can be invariantly dehned as a partial connection along the 
Hamiltonian vector held generated by the principal symbol (see [32], Section 3.1). 


Remark 4.1. Let M be a closed manifold and let A : C°°{M) C°°{M) be a 

classical pseudodiherential operator of order m. Some authors ([42]) like to identify 
C°°{M) with C'°°(M, H 1 / 2 ) by hrst introducing a smooth positive density p on M 
then identifying functions with half-densities by / -H- f^/p. Suppose this is the case 
then the sub-principal symbol of A locally dehned by ([42, Remark 2.1.10]) 


(4,3) 


Sub(4l) = a" + - . 


d‘^(j 


( 0 ) 

A 


+ 2 ' 


daf d{\ogp) 


dik 

is a homogeneous smooth function of degree m — 1 on T*M\0. We point out there 
is no essential diherence between (4.2) and (4.3) as it is always possible to endow 
M with a Riemannian metric g such that pg = p. To see this one can hrst endow 
M with an arbitrary Riemannian metric 'g, then get a positive smooth function k 
on M so that p = np-g, and hnally dehne g = where d is the dimension of M. 


Example 4.2. Let P = denote the connection Laplacian generated by the 
connection V on a vector bundle E of rank r. For hxed local frame for E\u 


X 

'( 0 ) 


1 




there exists a matrix oj = {oj 
Note in any local coordinates system 

Letting Vj = dj+hj one can easily get a'^p = g^’^{x)Q^k, CTp'* = —‘^'^g^’^bj^k+ig^^^fn^k- 
Considering = F^„ (see e.g. [5, Prop. 2.8], [ 6 , Sec. 2.5], [ 8 , Sec. 6 ]) we have 


of one-forms on U such that 0 

onU,P = - Ff^.Vfc) 

d) — C. I k^ln-rk 


Sub(P) = a^p^ -f 




( 0 ) 

p 


daP (9(log y/G) 


= a 


( 1 ) 


+ i 


dx^d^k 

dgp 


d^k 


dx^ 


dx’^ 


?. + i ■ g‘X 


-pn 

kn 


= + ig'Xnii.) + i' (-rLs”" - Tig"’) ■ + i ■ !f% 


■pn 
^ kn 


= -2ig^%Ak- 


Also, it is easy to check that hj = where denotes the transpose of oj. Thus 

in an invariant manner, Sub(P)(a:,.^) = —2ig{dx^,dx^)u'^{dj)^k = —‘^^gi^'^,0- 


Example 4.3. Let D = 7 V be the associated generalized Dirac operator of a Dirac 
bundle (P, 7 , V) of rank r. We adopt all the notations about V used in the previous 
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example. Then a® = 17 -^and consequently, 


Suh{D) = ^ ■ 

0Vg> i 87“> 

dilogVG) 

dx’^d^k 2 d^k 

dx^ 

= b - — — 
1 ^3 2 

d'y^ 1 u 

■ ■ 7 ^• FH 

dxk 2 


— b' — — 

70 ^ 2 

•(LtM-rLT”)- 

tOl H- 

'-1 

s ^ 

'I’^bk + bk'y^ 



where the third equality follows from (3.2) (see also the next subsection). On each 
inner product space (TxM,g^), we introduce two bilinear maps sending 

Xxj Yx G TxM respectively to ''y{Xx)uF{Yx) and {Xx)'^{Yx). Then it is easy to see 
that 


(4.4) 


Sub(Zi))(a:, ^) 


Tr( 4) + Tr(J4) 
2 


which means Sub(T))(a:, .^) actually is independent of 


With the sub-principal symbol concept available, we can express certain Wodzicki 
residues in a more invariant way. A key ingredient (see e.g. [13, 49]) is that if / is a 
smooth homogeneous function of degree 1 — d on then ■§^dS{^) = 0 for 

each k. For example, if A : E) —)■ E) is a classical pseudodifferential 

operator of order 1 — d, then it is easy to see that 



Tr(<77(l.O)o!SK) 



'n{Suh{A){x,0)dS{0, 


and consequently, the global density TeSx{A)dx is also of the form 


1 

( 2 ^ 


Tr(Sub(A)(x,0)t^*5(0)f^^, 


which is independent of the choice of local coordinates of M and local frames of E. 
For this reason, if A is a classical pseudodifferential operator of order 1 — d, then 
res (A) can be written as 

(4.5) res(A) = f Tr(Sub(A)), 

( 271 )“' Jt^m 

which is simply short for 


'M 


(2ir) 


'l£l=l 


Tr(Sub(A)(a:, .^))(iS'(.^) ]dx. 


Similarly, any integration over the unit cotangent bundle in the paper is always 
understood in this manner. 
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Proposition 4.4. Let A, B he classical pseudodifferential operators on sections of 
E such that the sum of the order of A and B is 1 — d. Then 

res(74i?) = [ Tr(Sub(A) ■ ■ Sub(i?)). 

( 27 r )'^ Jt*m 

Proof. It is straightforward to verify that (see also [10, (1.4)]) 

(4.6) Sub(y4il) = Sub(A) ■ a® + a® ■ Sub(i?) + — 


where 


Pa, A) = 


3af 



aA 

dik 

dx^ 

dx^ 

dik ' 

f} + {A\ 

A’}) = 

0, one c 


Tr(Sub(ylB) + Sub(aA)) ^ ^ 

which proves the proposition via (4.5) as res(74i?) = res(i?A). We are done. □ 


According to Proposition 2.8 and Proposition 4.4, one has 

(4.7) £/i{A,P) = j Tr(Sub(A) • • Sub(P'^^)). 

[Zti) Jt*m p ^ 

For any real number q, it is not hard to verify that (see also [10, (1.3)] and its proof 
therein) 

(4.8) Sub(P'') = q ■ ■ Sub(P). 

Thus by combining (4.7) ~ (4.8) and by considering dp, = (cip^)'^, we get 

Theorem 4.5. Let A be a classical pseudodifferential operator of order m and let 
P be a non-negative self-adjoint Laplacian on sections of E. Then 

(4.9) M(A,P) = -^/' Tr(Sub(A) + ^—■ Sub(P)). 

Jt*M ^ 


Theorem 1.1 is an immediate consequence of Theorem 4.5 as it suffices to first 
note from (2.24) that 


(4.10) 


^M,D) 


M(71,P2) + M(71P,P2) 


then apply Theorem 4.5 accordingly. 


Theorem 1.2 is also a consequence of Theorem 4.5, which can be seen as follows. 

Step 1: Let P be a smooth bundle endomorphism and let P be a self-adjoint 
Dirac type operator on sections of E. Obviously, it follows from Theorem 4.5 that 
^/i(P, Dy = 0. So according to (4.10), we get £^i{E., D) = gy Proposition 

2.8, we also have 

s^i{ED,Dy = ies{ED\D\-y = ies{E\D\-^D) = res(PP|P|-‘^) = s^i{DE,Dy. 
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Consequently, 

(4.11) 


A{F,D) = 


£/i{FD,D‘^) + £/i{DF,D‘^ 


Application of Theorem 4.5 results in 


^i{F,D) = 


(2vr) 


IT,*M 


Tt F- 


Sub(D) d • Sub(T)2) + Sub(D2) . a''^ 
2 


( 0 ) 


4 


Step 2: Recall D = yV + "0 is a self-adjoint Dirac type operator with potential ^ 
associated with the Dirac bundle (i?, 7 , V). In local coordinates we assume = 
dj + bj, where bj are smooth matrix-valued functions. Thus cr®(a;o,0 = i 7 '^( 3 ^o)'Ci 5 
cr}}^(a;o,0 = l^{xo)bj{xo) + ^jJ{xo). The compatible condition (3.2) between 7 and 
V gives = [ 7 ^, bj] — r^^ 7 "'. We further assume the local coordinate system is a 
Riemannian normal one centered at Xq G M. Hence 

(^)(^o) = [y\xo),b^{xQ)]. 

So by the definition of the sub-principal symbol of D, we get 

(4.12) Suh{D){xo,0 = 'y^ixo)bj{xo)+ij{xo)--{-^){xo) = C ^ 2 +V^)(^o)- 

Similarly, one can show that a^ 2 (a;o ,0 = + -‘^h){xo)^k- 

So by the dehnition of the sub-principal symbol of we get 

Snh{D^){xo, 0 = 1 ( 7 ^^ + '^1’" - ‘ 2 bk){xo)^k- 

Consequently, 

(4.13) (ag^ ■ Sub(D2))(xo,0 = -(7^7^'0 + 7'^'07'' - ‘27^h){xo)^j^k, 

(4.14) (Sub(D2) ■ ag^)(xo,0 = -(7^'07^ +-07^7^ - 2bk'y^){xo)^k^j. 


Step 3: By the above formulae we obtain 


[ Tr(F ■ Sub(D))(xo,0^^^(0 = Tr(F^^^^^±^ + F^l:){xo) ■ Vol(§"-i), 

7|q=i 2 

[ Tr(F ■ a® ■ Snh{D^)){xo,OdS{0 = Tr(FV' - ^ + 2^^){xo) • Vol(§"-'), 

Jg|=i « a 

[ Tr(F ■ Sub(D^) • a‘^^'>){xo,OdS{0 = Tr(FV' “ ^ + 2^^)(xo) ■ VoKS""'). 

Jg|=i « « 

Therefore, by the last formula in Step 1, 


which proves Theorem 1.2 as Vo^S"^"^) = 


i/j — {d — 2)ip 
4 

r(d/ 2 ) • 


) ■ Vol(§'^-^), 
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5. Dirac type operators of vanishing second coefeicient 


Let D = 7 V + ^ be a Dirac type operator of potential ^ associated with the Dirac 
bundle (-S', 7 , V). In this section we will no longer assume the self-adjointness of D. 
According to Remark 3.2, G C°°{M;End{E)) is well-defined. 

Theorem 5.1. Let D be a Dirac type operator. Then = 0 if and only if 

D is a generalized Dirac operator. 


By considering (1.16) and (1.17), we see that Theorem 5.1 implies Theorem 1.3 as 
a self-adjoint Dirac type operator is a Dirac type operator. But we should remark 
that Theorem 5.1 and Theorem 1.3 actually are equivalent as for any Dirac bundle 
one can always introduce a hermitian structure so that the associated generalized 
Dirac operator is self-adjoint on the Hilbert space defined by this hermitian structure 
and the metric measure on M (see e.g. [3]). 


Now we explain how to prove Theorem 5.1. Let D be a Dirac type operator of 
potential if associated with the Dirac bundle (77,7, V). According to Remark 3.2, 
J^(D, D^) = 0 if and only ii if = {d — 2)if. Thus to prove Theorem 5.1 it suffices to 
show that ii if = {d — 2)if, then there exists a connection V on T' compatible with 
7 such that D = 7 V. Note any connection on E must be of the form Vl = V -|- L, 
where 

L : C°^{M;TM) C'°°(M; End(E)) 

is a C°°(M)-linear map. It is easy to check that Vl is compatible with 7 if and only 
ii L{X) commutes with 7 (R) for all X,Y E C°°{M;TM). Suppose we do have such 
a map L such that D = 7 V l- Letting be a local orthonormal frame in TM, 

we have D = 7 V -|- 7 (efc)L(efc), which means if = j{ek)L{ek). Consequently, for any 
fixed i G {1 ,..., d}, 


-f{ei)if + if-f{ei 


= = _i(„ 


2 ' ' 2 
This means that globally L must be uniquely of the following form 

^{X)if + if^{X) 


(5.1) 


L{X) = 


for any X G C°°{M]TM). Thus to prove Theorem 5.1, with assuming if = {d — 2)if 
and defining V = V-l-L with L given by (5.1), we need only to show that 1) D = 7 V, 
and 2 ) V is compatible with 7 . 

The first property can be verified in the following way. Let {ek}f^i be a local 
orthonormal frame in TM, then 


7(efc)L(efc) = -7(6*,) 


'y{ek)if + if^ick) dif-if 


= V', 


2 2 
which gives D = 7 V + if = 7 V -|- 'y{ek)L{ek) = 7 V. 

To prove the second property we need to show that L{X) commutes with 7 (R) 
for all X,Y E C°°(M;TM). But considering L and 7 are C'°°(M)-linear maps, it 
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suffices to prove that Lx{Xx) commutes with '^xiXx) for ah Xx-,Yx G T^M at each 
X G M. This is guaranteed by the next proposition and thus a proof of Theorem 
5.1 is achieved. 


Proposition 5.2. Let (hh, 7) be a complex Cl{W^)-module. Let'll! G End(hT) be such 
that ijj = [d — 2)'ili, a'nd let L : —)■ End(hh) be the Imear map defined by 


UX) 




Then L{X) commutes with 7 (E) for all X,Y E 


The rest of this section is mainly devoted to proving Proposition 5.2. In the hrst 
subsection we begin with reviewing a few basic concepts about the representations of 
irreducible complex Clifford modules (see e.g. [45]), then give a proof of Proposition 
5.2. A few examples such as the massless Dirac operators ([ 8 , 9]) are studied in the 
second part. 


5.1. Clifford modules. Let C1(R'^) be the complex Clifford algebra generated by 
(•, •)) subject to the commutation relation X*Y-\-Y*X = —2{X, Y), where 
X,Y E (•, •) denotes the standard Euclidean metric on and * denotes the 
Clifford algebra operation. Any complex Cl(M'^)-module W with structure 7 (unital 
algebra morphism from C1(M'^) to End(hP)) studied in this paper is always assumed 
to be of dimension in N. A complex Cl(M'^)-module (hP, 7 ) is said to be irreducible if 
for any decomposition W = Wi © W 2 into subspaces invariant under 7 one has Wi = 
W or W 2 = W. It is known that a complex Cl(M'^)-module is irreducible if and only 
if it is of complex dimension 2 I- 2 J. Two complex Cl(M'^)-modules (lPi, 7 ) and ( 1 ^ 2 , 7 ) 
are said to be equivalent if there exists an invertible element k E Hom(lPi, IP 2 ) such 
that 7 = A* 7 , where n* is defined by sending q G End(lPi) to KgK,~^ G End(lp 2 )- Up 
to isomorphism there are exactly inequivalent irreducible complex C1(M'^)- 

modules. It is also known that any complex Cl(M'^)-module is a sum of irreducible 
ones. 

Let (VP, 7 ) be a complex Cl(M'^)-module and define -0 = Ylt=i for any 

ip E End(lP), where {ek\j.^i is an arbitrary orthonormal basis for (•, •)). Let 

Endfc(lP) := Span^{7(6*1 )7(ei2) •' -7(64) : 1 < A < *2 < • • • < 4 < d}, 

for fc > 0 and let Endo(lP) be the complex span of the identity map. A simple com¬ 
putation using the Clifford algebra relations shows that Endfc(lP) is an eigenspace for 
the map 'ip i-E ip with eigenvalue {—lfi{2k — d). It is clear that Clifford multiplication 
by the volume element 7 ( 61 ) 7 ( 62 ) • • - 7 ( 6 ^) defines a linear map from Endfc(fP) to 
Endrf_fc(fP). If the module is irreducible the subspaces Endfc(lP) generate End(lP) 
and therefore this gives a decomposition into eigenspaces. Hence, in the irreducible 
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case we have 


End(hh) = Endfc(hh) if d is even, 


k=0 


End(PE) = 0 EndkiW) if d is odd. 


k=0 

In the latter case we have used Endd-fc(hE) = Endfc(hE). This can be seen direclty 
because Clifford multiplication by the volume element commutes with the Clifford 
action. Therefore, by Schur’s lemma, it must be a multiple of the identity. As an 
immediate consequence of this discussion we get the following Proposition. 

Proposition 5.3. If (1T,7) is irreducible, then the d — 2 eigenspace of the map 
equals Endi(iy) = Spanc 7 (M'^). If furthermore d is odd, then 2 — d is not 
an eigenvalue of the map fj fj. 

Now we can start to prove Proposition 5.2. Let (IP, 7 ) be a (finite-dimensional) 
complex Cl(R‘^)-module and let (lPo, 7 ^°^) be an irreducible complex Cl(M'^)-module. 
It is easy to see that up to isomorphism (lPo, 7 ^°^) can represent all irreducible 
complex Cl(M'^)-modules if d is even, and so for (lTo,± 7 ®) if d is odd. We then 
have two cases to consider. 

Case 1: Suppose d is even. In this case there exists a (hnite-dimensional) complex 
vector space V such that (hP, 7) is equivalent to (IPo <8) E, 7), where the Clifford 
structure 7 on hPo < 8 ) E is applied only to the first tensor factor hPo. To be precise, 
7 (E) = 7 ^°^(E)( 8 )Idy for any Y G Without loss of generality we can assume that 
(VP, 7 ) = (hPo C) E, 7 ) since the general cases can be dealt with simply by studying 
isomorphism between Clifford modules. Then as an application of Proposition 5.3, 
we claim that the d — 2 eigenspace of the map ip equals 7 (°)(M'^) (g) End(E). To 
prove this claim, it suffices to hrst fix a linear basis {Ki\ for End(E), then express 
Ip with Ip = {d — 2)ip uniquely as ipf'^ ® Ki where ipf''^ G End(lPo), note 

® Ki = Ip = {d - 2)ip = {d- 2)ipi (g) Ki, 

and hnally apply Proposition 5.3 appropriately. This claim immediately implies that 

L{R<^) = {Idwo}(2)End(E), 

which obviously commutes with 7 (M'^) = 7 (°^(M'^) (g) {Idy}. 

Case 2: Suppose d is odd. In this case there exist (hnite-dimensional) complex 
vector space Ei, V 2 such that (hP, 7 ) is equivalent to ((lEo<g)Ei)©(lEo<g)E 2 ), 7 ), where 
the Clifford structure 7 on (hPo © Ei) © (hPo © V 2 ) applied only to hPo, is dehned by 

7(y) = (7'“>(y) ® Wh) © (- 7'“>(r) ® Wk,) 

for any Y G Similar to the discussions in the previous case, we assume without 
loss of generality that (hP, 7 ) = ((hPo © Ei) © (hPo © ¥ 2 ),^). Then as an application 
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of Proposition 5.3, we claim that the d — 2 eigenspace of the map ifj i/j equals 
(g) End(l/i)) © © End(V 2 )). To prove this claim, we hrst choose a 

linear basis for each Hom(V'Q,, Vjs), then express '0 with 'ip = (d — 2)'0 uniquely 

as 

(V-© 0 Af‘>) 0 0 Kf) + (i,r 0 Af 7 © (©'> 0 Af 7. 

where ^ End(lTo)- Here to be clear, W = {Wq © Vi) © {Wq © VP) is naturally 
identihed with (lTo®h 2 )©(hfo©Pi) by interchanging the positions, thus any element 
in W mapped via the second summand is indeed contained in W. Following this 
convention, it is straightforward to check that 

f = (7^ 0 Af 7 © (ir ® Af7 + (- V© 0 A<‘7 © (- vT 0 Af 7- 

Thus the claim is an immediate consequence of Proposition 5.3. This claim implies 
that 

= {{Idwo} ® End(Ei)) 0 {{IdwJ ® End(E 2 )), 

which commutes with 7 (E) for any Y G M'^. This hnishes the proof of Proposition 

5.2. 


5.2. Examples. In this part we study a few examples. Recall the diagonal values 
of the smooth kernel of De~^^ admit a uniform small-time asymptotic expansion 

CX) 

K{t,x,x,D,D‘^) ~ (f ^ 0+), 

k=0 


where ^ 

1 'ip — {d — 

(47r)'^/2 ^ 

If D is further self-adjoint, then the mollihed local counting function of D has a 
pointwise asymptotic expansion 




(x»/viXfi) ~ (a^ «3), 

k=0 


where 


Mi{D D'^) 

r(f) 


Example 5.4 (Generalized Dirac operators). Let D be a generalized Dirac operator 
associated with a Dirac bundle. Obviously, JYi{D,D^) = 0. To compare, it was 
shown in [3] that Tr(J^(D,D^)) = 0. 


Example 5.5 (Three-dimensional manifolds). Let D be a Dirac type operator of 
potential ip associated with a Dirac bundle of rank two over a three-dimensional 
closed Riemannian manifold M. We claim D^) = 0 if and only if Tr('0) = 0. 

To prove this claim, we note that if W is an irreducible complex Cl(R^)-module, then 
W is of complex dimension 2 and End(IE) = Endo(IE) © Endi(IE). Consequently, 
Endi(iy) is precisely the trace-free part of End(IE). 
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Example 5.6 (Parallelizable manifolds). In this example we construct generalized 
Dirac operators acting on E), where M is a closed parallelizable Riemannian 

manifold of dimension d, E = M x with rank r = 2 L 2 J. The Riemannian metric 
on M is denoted by g and the Riemannian connection on TM is denoted by V. 
Since M is parallelizable, there exist smooth real vector helds such that 

is an orthonormal basis for {T^M^gx) at each x G M. Let {Rk}f.^i be 
hxed complex matrices of size r x r with RjRk + RkRj = for all 1 < j, k < d. 

Then it is straightforward to verify that 'j : TM ^ End(E) dehned by 

(5.2) j{X) = g{X,Xk)Rk 

is an irreducible Cl(TM)-module structure. Dehne a flat connection V on E by 


01 


(5.3) 


Vx0 = Vx 


X 01 

X^r 


and dehne L : TM ^ End(E) by 

(5.4) UX) = 

Following Branson-Gilkey ([3, Lemma 1.3]), we claim V = V + L is a compatible 
connection on {E,'j). To this end we hrst let X,a E C°°{M-,TM), 0 G C°°{M; E). 
Obviously, a is of the form a = akXk for some ak G C°°{M). Thus 


Vx(7(«)0) = Vx(afcRfc0) = Vx{o'kRk4>) + L{X)'y{a)(j) 
= {Xak){Rkct>) + 7(a)Vx0 + L{X)-f{a)ct> 


and 

7(Vxa)0 + 7(«)Vx0 = 'i{{Xak)Xk + akV xXk)(j) + 7(q!)Vx0 + 'y{a)L{X)(j) 

= {Xak){Rk(l>) + afc7(VxXfc)0 + 7(a)Vx0 + 7(a)^(^)0- 


Hence to prove the compatible condition (3.2) currently for a G C°°{M]TM) it 
suffices to show [L(X), 7 (a)] = g{a, Xi:)'j(Vx^k), which is obviously equivalent to 

(5.5) [L(X,),7(X,)] = g{Xj,XMXxM (1 < t,j < d). 

Letting Th = —L^^ denote the Christoffel symbols, that is, XxiXj = T^jXk, we have 

[i(X),7(A)l = \(Rk7Xx.Xt)R, - RjRk'tXx.Xk)) 

= \ {RX7kRnRj - RjRtrttRn) 

= 1 ( - - RkTlR,R,^ - RiRtVlkK) 

= \i-RXl + niRn) = nfin 

= iXx.Xi) = g(.Xj,Xth(.Vx.Xk), 
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which proves (5.5). Generally, if a G Cl(TM)) is of the form a = * 

■ • ■ * where G C°°{M;TM) {1 < k < n), then according to (3.1), 


Vx(7(a)0) = Vx(7(«^^^) 7(«^^^) • • •7(a^"^)0 ) 

= 7 (Vx«’'^^) 7(0*'^^) • • •7(a^"'^)0 +7(a*'^^)Vx 7(0*'^^) • • • 7(a*^'^^)0 


n 

= ^7(0;^^^) • • • 7(a*-^“^^)7(Vx«*'^^)7(«^^’'’^^) • • • 7(a*'"'^)(/> + 7(a)Vx0 

k=l 

= 7(Vxa)0 + 7(«) Vx0, 


which snffices to claim (3.2) for arbitrary a G G°°(M; Cl(TM)) by linearity. Thus 
(E, 7 , V) is a Dirac bundle and 


(5.6) 


D = 7 V = 7 V + 


RiRjiiy Xi^j 


is the associated generalized Dirac operator. By Theorem 5.1, we get Mi{D,D‘^) = 
0. It is well-known ([3]) that one can always furnish (T^, 7 , V) with a smooth her- 
mitian hber metric so that D is self-adjoint on L‘^{M] E) dehned by this fiber metric 
and the metric measure on M. So in this context ^i{D) = 0. Next we derive a 
local expression of D. To this end we first dehne a flat connection V on End(£') by 


(5.7) 



then introduce 

(5.8) T(X, F, Z, W) = ^{XMY) {^{VzW) - Vz{l{W))), 

where X, Y, Z,W G C°°{M]TM). It is easy to verify that T is G°°(M)-linear in all 
of the four variables. This implies Tr 24 (Tri 3 (T)) is a smooth endomorphism of E, 
where Tri 3 (\h) denotes the trace of T with respect to the hrst and third variables, 
and Tr 2 , 4 (T) is understood in a similar way. Thus globally we have 


Tr 2 , 4 (Tri, 3 (T)) = ^{Xi, Xj, X^, Xj) 


RiRjiiy XiXj), 


which means the generalized Dirac operator D can also be written as 

Tr 2 , 4 (Tri, 3 (T)) 


(5.9) 


D = 7 V -I- 


Hence in a local coordinate system [x ^,..., x'^), D is of the form 
(5.10) D = 
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As an example, we assume that M is a three-dimensional oriented closed Riemannian 
manifold. By Steenrod’s theorem M is parallelizable and we can set 


Ri = 


0 -i 
-i 0 


R, = 


0 -1 

1 0 


R.. = 


-i 0 
0 i 


Then it is straightforward to verify that D is nothing but the massless Dirac operator 
W locally defined by (A.3) in [8]. Note Tte{^i{D)) = 0 was established in [8]. To 
compare, we have shown that = 0. Moreover, it was proved in [8] that the 

sub-principal symbol of D is proportional to the identity matrix at each point of the 
manifold M, where the corresponding bundle trivialization {si, S 2 } of E = M x 
is chosen to be 


Si = 


S2 = 


In the following we will give a short proof of this fact. According to Example 4.3 
and our previous discussions, it is easy to see that oj = L, and consequently, 

Sub(£>) = l(7(V)i(Y) + L(V,)7(V,)) 

= g{RjRkf^^kRn + Rkf^'jkRnRj) ■, 

where we recall the Christoffel symbols are given by Vx^Xk = with 

T”^ = —T^^. Note Sub(D) is a sum of 27 = 3 x 3 x 3 items. Next we no longer use 


Einstein’s sum convention and decompose Sub(D) into four parts: 

Sub(D) = ^x+ *+ A+ Y ^ = I + n + III + IV, 

k=n k^nj=n k^n,j^n^k=j k^nj^n^k^j 

where x is short for ^{RjRk f^'jkRn + Rkf^ ^i^RnRj)- One can easily deduce from the 
properties T"^ = —T^^, RjRk + RkRj = —‘25jk that I = II = III = 0 and 

IV = l{RjRkT],R,, + RkT],Rr,R,) 


{j,k,n)eiT3 

1 
4 


T],R^RkRn 


{j,k,n)&iT 3 


1 

2 


Y T^^RjRkRn 


{j,k,n)&iT3 is even 
r3 I ri I r2 

i 12 W j- 9 C[ W j- q 
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where tts denotes the symmetric group of degree 3, and the last equality follows from 
R 1 R 2 R 3 = R 2 R 3 R 1 = R 3 R 1 R 2 = — 1. To conclude we get 


(5.11) 


Sub(D) = 


n2 + r 


1 

23 


p2 

31 


which agrees with [8, Lemma 6.1]. Furthermore, it is characterized in [8] that given 
a self-adjoint Dirac type operator D on sections of M x C^, it is a massless one if and 
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only if = 0 and Sub(-D) is pointwise proportional to the identity matrix, 

where the corresponding bundle trivialization {si, S 2 } oi M x is chosen as before. 
We can also give a proof of this fact. Based on the previous discussions, it suffices 
to prove the necessary part. So let D be a self-adjoint Dirac type operator such that 
Tr(..^i(D)) = 0 and Sub(D) is pointwise proportional to the identity matrix. Since 
M is a parallelizable three-dimensional manifold, it is easy to see that there exists a 
unique massless Dirac operator 7 V sharing the same principal symbol with D and 
a unique bundle endomorphism ^|J such that D = 7 V -|- ip- According to Example 
5.5, Tr('0) = 0. Note also 

Sub(D) = Sub( 7 V) -I- ip, 

which implies that ip is pointwise proportional to the identity matrix. Thus we must 
have '0 = 0, and consequently, D = 7 V is indeed a massless Dirac operator. 
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